Dilepton production by bremsstrahlung of 
meson fields in nuclear collisions 



X 



I.N. Mishustin'^^ L.M. Satarov'^, H. Stocker^ and W. Greiner'^ 

The Kurchatov Institute, 123182 Moscow, Russia 
^ The Niels Bohr Institute, DK-2100 Copenhagen 0, Denmark 
Institut fiir Theoretische Physik, J. W. Goethe Universitdt, 
D-6OO54 Frankfurt am Main, Germany 

Abstract 



We study the bremsstrahlung of vh'tual cj-mesons due to the collective deceleration of 
nuclei at the initial stage of an ultrarelativistic heavy-ion collision. It is shown that elec- 
tromagnetic decays of these mesons may give an important contribution to the observed 
\ yields of dileptons. Mass spectra of e'^e~ and fi^ pairs produced in central Au+Au 

' collisions are calculated under some simplifying assumptions on the space-time variation 

of the baryonic current in a nuclear collision process. Comparison with the CERES data 
for 160 AGev Pb+Au collisions shows that the proposed mechanism gives a noticeable 
^ ' fraction of the observed e~^e~ pairs in the intermediate region of invariant masses. Sensi- 

'nI" . tivity of the dilepton yield to the in-medium modification of masses and widths of vector 

^ . mesons is demonstrated. 

a^ 
o 

' According to the relativistic mean-field model [||], strong time-dependent meson fields are 

r-| ■ generated in the course of a relativistic heavy-ion collision. Using the approach developed in 

p papers on the pion 0] and photon bremsstrahlung we suggested recently IQ a new mechanism 

' of particle production by the collective bremsstrahlung and decay of classical meson fields in 



relativistic heavy-ion collisions. The comparison with the observed data on pion multiplicity 



■ shows that this mechanism may be important in central collisions of most heavy nuclei 

already at the SPS bombarding energy. 

Within this mechanism the production of some particle (s) i is considered as a two-step 

process Ap + At ^ uj* ^ i + X . Here Ap{At) stands for the projectile (target) nucleus and uj* 

is a virtual vector meson Q. The first step in the above reaction corresponds to the virtual 

bremsstrahlung process leading to the creation of an off-mass-shell vector meson. The second 

step is the superposition of all channels of the virtual meson decay with the particle i in the final 

state. Below we consider the production of virtual mesons in the coherent process caused by 

the collective deceleration of the projectile and target nuclei at the initial stage of the reaction. 
^Here and below the virtual particles {u and 7) are marked by a superscript 
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Preliminary results concerning the contribution of the above mechanism to the production of 
pions, real cu-mesons and dileptons are published in Ref. 

Below we focus mainly on the dilepton emission in central collisions of ultrarelativistic 
nuclei. The analysis of the dilepton production is interesting at least by two reasons. First, 
dileptons are highly penetrating particles and carry practically an undistorted information 
about their creation points. Second, a strong enhancement of the dilepton yield was observed 
recently in central 200 AGeV S+Au 0, S+W and 160 AGeV Pb+Au [|] collisions. The 
analysis shows [9-11] that this enhancement can only partially be explained by the conventional 
mechanism of binary hadron collisions, e.g. by tttt ^ p ^ e+e^ processes. According to Refs. P, 



10] the experimental data can be reproduced assuming a strong reduction of the p meson mass 



in dense and hot nuclear matter. On the other hand, as argued in Ref. |TT]], the in-medium effect 
was probably overestimated in these calculations. Below we show that the enhanced dilepton 
yield may be explained, at least partly, by the contribution of the collective bremsstrahlung 
mechanism. 

1. Formulation of the model 

By analogy to the Walecka model we introduce the vector meson field uj^{x) coupled to the 
4-current J^{x) of baryons participating in a heavy-ion collision at a given impact parameter b. 
The equation of motion determining the space-time behaviour of co''^(x) can be written as 
{c = n = l) 

{d''d, + ml)LU^{x)=gvnx), (1) 

where gy is the uNN coupling constant and m^^ ~ 783 MeV is the omega meson mass. In the 
mean-field approximation is considered as a purely classical field. From Eq. (|l]) one can see 
that excitation of propagating waves in the vacuum (the bremsstrahlung process) is possible if 
the Fourier transformed baryonic current 



JA'(p) = / d^xJ^{x)e'P'= (2) 
is nonzero in the time-like region p'^ ~ 
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In the following we study the bremsstrahlung process in the lowest order approximation 
neglecting the back reaction and reabsorption of the emitted vector mesons, i.e. treating J'^ as 
an external current. From Eq. (0) one can calculate the energy flux of the vector field at a large 
distance from the collision region . Then this flux is expressed in terms of the distribution of 
field quanta, i.e. cj-mesons. This leads to the following formulae for the momentum distribution 



of real cu-mesons emitted in a heavy-ion collision [12 



d^A^ 

E^-^ = S{E^,p), (3) 
dp 



where E^ = Jm'^ + p'^ and 



sip) = -^smp)j'(p)\ (4) 

is a source function. In our model the latter is fully determined by the collective motion of the 
projectile and target nucleons. 

To take into account the off-mass-shell effects we characterize virtual u mesons by the mass 
= \fp^ and total width V^* . The spectral function of virtual oj mesons may be written as 



To calculate the distribution of virtual mesons in their 4-momenta p we use the formulae 

^ = p{M^)S{p) . (6) 
d p 

In the limit F^^* — > one can replace p{Muj) by 2(5(M2 — 'tiJ) • In this case Eq. (||) becomes 
equivalent to the formulae (|^) for the spectrum of the on-mass-shell vector mesons. Below we 
study also how the dilepton production is changed when the pole position in the vector meson 
propagator is shifted due to the medium effects. 

We consider the following channels of the virtual a; decay, most important at invariant 
masses M^<m^: i = e"'"e~, 7r°7, 7r°e~'"e~, 7r°/i+/i~, tt+tt", n^n^n^ . The total width F^^* 

is decomposed into the sum of partial decay widths F^*^j : 
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The distribution over the total 4-momentum of particles in a given decay channel can be written 



as 



-4 - ^u,*~*i-74 ' l°J 

dp dp 

where B^*^^ = r^*^^/ri^* is the branching ratio of the i-th decay channel. The latter is a 
function of the total invariant mass of the decay particles M = ^Jp^ = iW^. 

To calculate the 4-current J'^{p) determining the source function S{p) we adopt the simple 
picture of a high-energy heavy-ion collision suggested in Ref. [Q. We consider collisions of 
identical nuclei {Ap = At = A) at zero impact parameter. In the equal velocity frame the 
projectile and target nuclei initially move towards each other with velocities ±fo or rapidities 
±?/0) where Vq = tanhyo = (1 — 4m^/s)^/^ and ^/s is the cm. bombarding energy per nucleon. 
In the "frozen density" approximation |Q the internal compression and transverse motion of 
nuclear matter are disregarded at the early (interpenetration) stage of the reaction. Within 
this approximation the colliding nuclei move as a whole along the beam axis with instantaneous 
velocities Zp = —Zt = z{t) . The projectile velocity z{t) is a decreasing function of time, which 
we parametrize in the form H] 



where r is the effective deceleration time and f/ is the final velocity of nuclei (at t — > +00) . 

In this approximation the Fourier transform of the baryon current J'^{p) is totally deter- 
mined by the projectile trajectory z{t) 

00 

J0(p) = ^ J3(p) = 2A f dte'P^^ cos \p,z{t)] F (Jp^ + p^ . [I - z^{t)]) , (10) 

Po J 

—00 

where and ^^^^5 respectively, the longitudinal and transverse components of the three- 
momentum p, F{q) is the density form factor of the initial nuclei 

1 ■ 4 °° 

F(q) = — dc'rp{r)e~^^'^ = — / rdrp(r) sin gr . (11) 

/\ J /\Q J 



The time integrals in Eq. ([T0|) were calculated numerically assuming the Woods-Saxon distri- 
bution of the nuclear density p(r) . 
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In this work we choose the same couphng constant, (y'i/=13.78 and stopping parameters r, Vf 
as in Ref. [Q. In particular, it is assumed that r equals one half of the nuclear passage time 

it 

r=— , 12 
smh I/O 

where R is the geometrical radius of initial nuclei, R = tqA^^^ with tq = 1.12 fm. Instead of f/ 
it is more convenient to introduce the cm. rapidity loss 6y defined by the relation: 

Vf = tanh{yo~ Sy) . (13) 



For central Au+Au collisions we assume the energy-independent value |T3| 6y = 2A for y/s > 10 
GeV and the full stopping {6y = y^) for lower bombarding energies. 

2. Different channels of virtual meson decay and dilepton 
production 



Similarly to Ref. |]T4l we assume that the "direct" decays of uj mesons into dileptons, 
uj* , proceed via the intermediate emission and decay of virtual photons 7*. In the 

following the explicit formulae are presented for the e"*"e~ production only. The corresponding 
expressions for dimuons are obtained by replacing the lepton masses and decay widths. The 
matrix element of the process u* e'^e^ is proportional to the u meson polarization vector 

the lepton current v+'^^U- and the photon propagator /c^^, where A; = p+ +p_ is the total 



4-momentum of the lepton pair. The calculation of the partial decay width gives the result [14 



r,*_.,,(M)cxM-^r,*^,,. (14) 

Here M = y/k^ is the dilepton invariant mass (in the direct channel M = M^J) and r^*^gg is 
the partial width of a virtual photon: 

r,*^ee = ^ (1 - y) - 2^e) , (15) 

where a = e^/^c, 9(x) = |(1 + signx), me is the electron mass and 



The proportionality coefficient in Eq. ( ]T^ ) is determined from the condition B^*^^i,{m^) = B^e 
where B^e = 7.1 ■ 10~^ is the observed probabihty of the — ee decay [P^. 

We take into account also the three-particle, "Dalitz" decays uj* n^e^e^ . At fixed values 
of M and the components of the total dilepton 4-momentum k in the u rest frame can be 
found by using the expressions 
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+ (17) 

where m^r = 0.14 GeV is the pion mass. The corresponding partial width can be calculated 
assuming that the Dalitz decay is the two-step process u* —>■ 717* Tree. Generalizing the 



results of [110] to the case of virtual cj's we obtain the following expression for the differential 



width of the Dalitz decay 



The UJ* — > irj* decay width is proportional to the electromagnetic form factor squared, F^^: 

r^*^,^* cx F^JkfQ{M^ -M-m^), (19) 

where \k\ is determined from Eq. (|13). The coefficient of proportionality may be found by 
considering the limiting case m^, M ^ 0, when the left hand side of Eq. ([T^) coincides 

with the observed width of the — > 7r7 decay. As shown in |M the experimental data for 



are well reproduced within the vector meson dominance model [T^l . Assuming that this model 
is valid also for decays of virtual cj's we have 

ml (ml + r^) 

where rrip and Tp are the mass and total width of the p meson. Unless otherwise stated, Eq. (^) 



is used with the parametrization Tp = Tp{M) suggested in Ref. ||Tj] and the free p meson mass 
{rrip = rripo ~ 0.77 GeV). 

To calculate the total width of virtual cu's one should also know the partial widths of non- 
electromagnetic decay channels. In the considered region of masses < 1 GeV we take into 
account the decays with two and three pions in the final state. Assuming that the uj* —>■ 27r 



6 



matrix element is proportional to the product of the uj meson polarization vector and the 
difference of the pion 4-momenta, we get 

r.*^2.(Mj oc hC ■ [Ml - 4m^)'^' e(M^ - 2m.) . (21) 



The proportionality coefficient is taken from the condition -B(j*^27r(^w) = -6277 = 0.022 [p!5| 



The 10* Stt partial width is calculated assuming that it is proportional to the 3-pion phase 



space volume allowed by the kinematics |12|: 



r^*-.3.(M.) « 'f3.(M,) ■ 0(M^ - 3m.) , (22) 



with the coefficient determined from the relation i?^*_^3.(mj^) = i?3. = 0.89 |TI 

The mass distribution of e+e~ pairs produced by the bremsstrahlung mechanism can be 
written as a convolution of the virtual lo meson spectrum and the differential branching of the 
LO* — > eeX decay {X denotes any particle(s) emitted together with the lepton pair): 



dM 7 ^ ■ dM ■ 

Taking into account only the direct and Dalitz decays, the dilepton mass distribution can be 
represented as: 

00 

dM- = ^-*-dM-+ I — • ('^^ 

The first (direct) term is obtained by using the relation dB^^*^^JdM = B^*^^^S{M — M^) and 
performing the explicit integration over M^. 

To compare the model predictions with experimental data one should take into account 
the various acceptance cuts used in different experiments. This severely complicates the cal- 
culations: in general one should know the differential branching d^Buj*^eex/d^p^d^p_ which 
describes the probability of the u meson decay into the lepton pair with the positron and 
electron 3-momenta and p_, respectively. To calculate the acceptance weighted mass dis- 
tribution one should replace di?^*^eex/dM in Eq. (|23|) by 

■ - ^ d'p+d'p. A ■ y,*r''' 6{M - v^) . (25) 



dM 7 ^ d^p+d-'p^ 
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The weight function A equals one (zero) if are inside (outside) the kinematical volume 
covered in a given experiment. In numerical calculations presented in the next section we use 
for e~^e~ pairs the acceptance cuts of the CERES experiment [§] 

Pt± > 0.175 GeV/c, 2.1 < r/± < 2.65, O^e > 0.035 . (26) 

Here pt± and ri± are the transverse momenta and pseudorapidities of leptons, 6ee is their relative 
emission angle in the lab frame. 

At given M and the components of the vectors are fixed by the angular variables 
Q = {6, 0) and Q = {6, 0), where Q denotes spherical angles of k with respect to the u meson 
3-momentum (in its rest frame) and Cl stands for the positron emission angles with respect to 
k (in the pair rest frame). By using these variables one can rewrite Eq. (^) as follows 

dM "dM J J dQdQ ' ^ ^ 

where dW^* ^^ex / ^^^^ is the angular distribution of the u* ^ eeX decay normalized to unity. 

In our case the usual procedure of calculating the direct and Dalitz parts of dilepton distri- 
butions by simple averaging over all polarizations of decaying u mesons is not correct. Indeed, 
as seen from Eq. (|T]), the polarization vector of a virtual vector meson is proportional to 
Jfj,{p) where p is the meson 4-momentum. In the uj meson rest frame = (0,^), where ^ is a 
vector parallel to the direction of p. Proceeding from the u* ee matrix element (see above) 
we get the following relation for the direct part of the angular distribution 

= C^. (l - cos^ ~e) , (28) 

where Cdir is the normalization constant. The anisotropy of the lepton angular distribution in 
the rest frame of the u meson is a consequence of its polarization. 



The Dalitz part of the dilepton distribution is calculated assuming ||TJ, |T6[ that the u* een 
matrix element is proportional to ^^uasi^V^^" ^+1^'^-) ■ The direct calculation gives 



dW„* 



u) ^ee-K ^ •2/1 

= CHa.i sm 



/?2sin2^sin2(0-0) , (29) 



dmn 

where Cdai is found from the normalization condition. The averaging over the uj meson polar- 
izations is equivalent to the averaging over VL. As a result we obtain the distribution over VL 



obtained earlier in Ref. 18 . 



3. Results 



Let us now discuss the results of numerical calculations obtained within the model described 
above. One should bear in mind that the model assumes a rather simplified space — time evo- 
lution, in particular the collective projectile-target deceleration (see the discussion in Ref. 0). 
Therefore, the model in its present form can be used for a qualitative analysis only. 

Fig. 1 shows the dilepton mass spectrum in central (6 = 0) 160 AGeV Au+Au collisions. 
One can see that the mass distribution of dileptons produced by the virtual uo decays is quite 
different from that predicted by the conventional hadronic sources. In particular, the low and 
intermediate mass region is strongly enhanced. This is explained by the copious production 
of "soft" virtual cu's by the bremsstrahlung mechanism. The contribution of direct uo decays 
is peaked at very small invariant masses as well as at the pole position of the uj propagator, 
M = m^. The Dalitz contribution is most important in the intermediate region of dilepton 
masses, 0.2 GeV <M< 0.6 GeV. A similar behaviour is predicted for the dimuon spectrum 0, 
Fig. 2. The main difference here is the much higher mass threshold at M = 2m^. 

In Fig. 3 we compare the model predictions for the same reaction, but at different bom- 
barding energies, ^/s = 17 A3 AGeV (SPS) and 200 AGeV (RHIC). We have also calculated the 
dilepton spectra at the LHC energy y/s=6.3 ATeV but the corresponding results practically co- 
incide with the model prediction for the RHIC energy. Such a behavior follows from the energy 
independence of the stopping parameter 6y, assumed at high i/i (see Sect. 1). As a conse- 
quence, the phase-space distribution of primordial lu mesons, produced by the bremsstrahlung 
mechanism, saturates with raising bombarding energy [Q. 

The experimental acceptance cuts and a poor mass resolution distort significantly the dilep- 
ton mass distributions as compared to those presented in Figs. 1-3. In Fig. 4 we show the dilep- 
ton mass spectrum for central 160 AGeV Au+Au collisions. The CERES acceptance cuts and 
mass resolution are included in this calculation. The double differential spectrum (PNee/ dMdrj 
is obtained by dividing the acceptance weighted mass distribution, Eq. (plD , by the width of 

the CERES pseudorapidity window. In calculating the acceptance weight A entering Eq. (^) 

^ Since the branching ratio of the direct decay uj is not known experimentally ||l5|| we assume that 

it is equal to B{uj e+e"). 
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we have neglected the transverse momenta of primordial uo mesons (see Ref. [Q). At the same 
figure we show separately the contributions of direct and Dalitz decays of uo mesons. Note that 
the original spectrum (without mass resolution corrections) of e"*"e~ pairs has a strong peak 
at M ~ TJiu)- This peak originates from the direct uj meson decays. The step-like behaviour 
of the direct contribution at M ~ 0.35 GeV appears due to the CERES cut at small trans- 
verse momenta, pt± > Pmin=0.175 GeV/c. Indeed, in the limit = the minimal invariant 
mass of "direct" pairs is 2 ml + p^j^ ^ 0.35 GeV. Taking into account nonzero components 
of prp^ will result in a certain smoothening of the above-mentioned jump in the dilepton mass 
distribution. 

As one can see from Fig. 4, the bremsstrahlung mechanism gives a significant contribution 
to the dilepton production in the intermediate mass region. Hopefully, this contribution can 
be observed experimentally by using a characteristic angular distribution in the dilepton mo- 
mentum predicted by the model 0, ^. However, the bremsstrahlung contribution alone is not 
sufficient to explain the dilepton yield observed in central 160 AGeV Pb+Au collisions. In the 
most interesting region of masses M ~ 0.4 — 0.6 GeV the data are underestimated by a factor 
of about three. 

A special calculation showed that this discrepancy can not be removed by taking into account 
the excitation and decays of virtual p mesons disregarded in the present calculation. The 
contribution of the p meson bremsstrahlung in symmetric heavy-ion collisions is proportional 
to the isospin asymmetry factor x = (1/2 — Z/A^ where Z and A are the charge and mass 
numbers of the colliding nuclei. Since x~10~^, the role of the p meson bremsstrahlung is 
relatively small even for heaviest nuclei. 

Of course, in addition to the collective bremsstrahlung mechanism the usual, incoherent, 
sources of dilepton production (e.g. the vrvr ee and vr — > ee7 processes) also give a noticeable 
contributions. According to Refs. [9-11] the dynamical models incorporating only these inco- 
herent hadronic sources may easily explain the low mass dilepton yields in S+Au and Pb+Au 
collisions at the SPS energies. On the other hand, these models significantly underestimate 
the observed data in the intermediate mass region. As argued in Refs. [|, |T0| the agreement 
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with experimental data can be achieved if one assumes a strong reduction of the vector meson 
masses in dense nuclear matter. 

To check the sensitivity of our model to the in-medium modification of the vector mesons, 
in Fig. 5 we compare the dilepton mass distributions calculated for different values of the p 
and uo masses. To diminish the number of model parameters we take the same mass reduction 
factor for p and uj mesons: mp/nipQ = m^/m^Q = R^. One can see that the dilepton mass 
distribution is rather sensitive to Rm- A relatively small, 20 %, reduction of the meson masses 
raises the dilepton yield at M ~ 0.5 GeV by a factor of about two. On the other hand, the 
model calculation with fixed Rm predicts too high peaks in mass distributions. One should 
bear in mind that these calculations provide only a rough estimate of the possible effect since 
in an actual nuclear collision mass shifts should be time (and space) dependent. Therefore, the 
observed distribution will be a superposition of contributions with different R^ ■ As a result, 
the peak of direct dileptons will be less pronounced in this distribution. 



As calculations of Ref. [1^ show, the cu meson width may be significantly increased in 
dense baryonic matter due to the mixing of a and uj mesons At baryonic densities of 
about two times the normal nuclear density, the width may increase by a factor of about 5 as 
compared to its vacuum value. To estimate the effects of the in-medium broadening of virtual 
UJ mesons, we have performed the calculation where the partial decay widths were scaled by the 
same amplification factor A, independent of M^. The calculation shows that this broadening 
influences mainly the direct component of the dilepton yield. As seen from Fig. 6, the yield in 
the intermediate mass region may be strongly enhanced at A>5. 

4. Summary 

In conclusion, we have shown that the collective bremsstrahlung of the vector meson field can 
provide an important source of dilepton production in high energy heavy-ion collisions. This 
mechanism may be responsible, at least partly, for the enhanced yield of dileptons observed 
in central nuclear collisions at the SPS bombarding energies. It has been demonstrated that 
the coherent dilepton production is sensitive to the in-medium modification of vector meson 
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masses and widths. Obviously, these effects should be studied in more details in microscopic 
models. 

In future studies one should implement a more realistic dynamical picture of a heavy-ion 
collision by using either the fluid-dynamical or kinetic approaches. In this way one can take 
into account the flow and compression effects disregarded in present model. Also the formalism 
should be generalized to study bremsstrahlung effects in a situation when masses and widths 
of vector mesons are space and time dependent. To extend the calculations to collider energies, 

> 200 GeV, the model should be reformulated at the quark-gluon level. 
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Figure captions 

Fig. 1. Mass spectrum of e^e~ pairs produced by tlie collective bremsstralilung mechanism 
in central 160 AGeV Au+Au collisions. Contributions of the direct (cj* — > ee) and 
Dalitz {ijj* — » Tree) decay channels are shown by dotted and dashed line respectively. 

Fig. 2. The same as in Fig. 1, but for the spectrum of fi'^fi" pairs. 

Fig. 3. Comparison of e^e~ yields in central Au+Au collisions at the SPS and RHIC bom- 
barding energies. 

Fig. 4. Mass spectrum of e^e^ pairs produced by the collective bremsstrahlung mechanism 
in central 160 AGeV Au+Au collisions. The grey histogram (solid line) shows the 
model results with (without) inclusion of experimental mass resolution. The dotted 
(dashed) line shows the contribution of direct (Dalitz) u meson decays. Preliminary 
experimental data for central Pb+Au collisions are shown by solid circles. 

Fig. 5. Comparison of e^e~ spectra in central 160 AGeV Au+Au collision for different 
values of the meson mass reduction factor Rm- Preliminary experimental data for 
central Pb+Au collisions are taken from Ref. p. 

Fig. 6. The same as in Fig. 5 but for different choices of the virtual u width in units of the 
vacuum width Fq = 8.4 MeV. 
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